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ABSTRACT 

Using high signal-to- noise ratio VLT/F0RS2 long-slit spectroscopy, we have stud- 
ied the properties of the central stellar populations and dynamics of a sample of SO 
galaxies in the Fornax Cluster. The central absorption-line indices in these galax- 
ies correlate well with the central velocity dispersions (ctq) in accordance with what 
previous studies found for elliptical galaxies. However, contrary to what it is usually 
assumed for cluster ellipticals, the observed correlations seem to be driven by system- 
atic age and a-element abundance variations, and not changes in overall metallicity. 
We also found that the observed scatter in the index-tJo relations can be partially 
explained by the rotationally-supported nature of these systems. Indeed, even tighter 
correlations exist between the line indices and the maximum circular velocity of the 
galaxies. This study suggests that the dynamical mass is the physical property driving 
these correlations, and for SO galaxies such masses have to be estimated assuming a 
large degree of rotational support. The observed trends imply that the most massive 
SOs have the shortest star-formation timescales and the oldest stellar populations. 

Key vifords: galaxies: elliptical and lenticular - galaxies: stellar populations 



This is the third paper of a series designed to study SO 
galaxies in the Fornax Cluster by using optical, long-slit 
spectroscopy (VLT/F0RS2) and archive optical and near- 
infrared imaging. In Bedregal et al. (2006a, hereafter Pa- 
per I) we described the sample and studied the stellar kine- 
matics along the apparent semimajor axes of these objects. 
In Bedregal et al. (2006b, hereafter Paper II) we use the cir- 
cular rotational velocities obtained in Paper I to study the 
TuUy-Fisher relation (TuUy & Fisher 1977) of these galaxies 
as part of a larger compiled sample of ~ 60 local SOs. In 
this paper we study the sample of Fornax SOs and try to 
establish links between the properties of their central stellar 
populations (ages, chemical abundances) and their global 
properties such as mass and dynamics. The ultimate goal 
is to find the main physical drivers governing the formation 
and evolution of SOs. 

Many previous works on stellar populations of early- 
type galaxies have concentrated on the relations between 
velocity dispersion and Mg indices (e.g. Terlevich et al. 
1981; Gorgas et al. 1990; Guzman et al. 1992; Bender et 
al. 1993; Jorgensen et al. 1996; Bender et al. 1998; Bernardi 
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et al. 1998; Colless et al. 1999; Jorgensen 1999; Kuntschner 
2000; Poggianti et al. 2001; Mehlert et al. 2003; Thomas 
et al. 2004), which are usually interpreted as correlations 
between mass and metallicity. However, there is still con- 
troversy about the role of relative ages and abundances of 
different elements in the observed trends (Jorgensen 1999; 
Trager et al. 1998, 2000; Kuntschner et al. 2001; Poggianti et 
al. 2001; Mehlert et al. 2003; Thomas et al. 2005) with some 
recent results pointing to an important dependence of these 
relations on both age and a-element abundances (Gallazzi 
et al. 2006; Sanchez-Blazquez et al. 2006). 

The vast majority of these studies bundle together el- 
lipticals and SOs treating them as a single population of 
'early-type' galaxies. Here we concentrate on SOs as a dis- 
tinct class. Their formation history could be very different 
from that of ellipticals even though some of their properties 
may appear similar (see, e.g., Aragon-Salamanca, Bedregal 
and Merrifield 2006). By using a full set of 10 absorption 
line indices and careful comparison with the results from 
stellar population synthesis models we will try to break the 
degeneracies between age, overall metallicity and a-element 
abundances. 

The remainder of the paper is laid out as follows. In 
Section [2l we describe the different steps followed to obtain 
accurate Lick absorption line indices. In Section [3] we make 
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some tests in order to compare Lick index measurements at 
different spectral resolutions. In Section |4] we make consis- 
tency checks between the data and the Bruzual & Chariot 
(2003) models. Section [5] presents the main results and then 
discusses their implications. Finally in Section |6] our conclu- 
sions are summarised. 



2 THE DATA: LINE INDEX MEASUREMENT 

In this section, the main steps followed to calculate Lick/IDS 
line-strength indices (Burstein et al. 1984, 1986; Worthey et 
al. 1994; Worthey & Ottaviani 1997) using long-slit spectra 
are described. This will allow us to study scaling relations 
such as the one between Mg2 and central velocity disper- 
sion (e.g. Burstein et al. 1988; Guzman et al. 1992; Bender, 
Burstein & Faber 1993) and to compare the results to pre- 
dictions from simple stellar population models of Bruzual & 
Chariot (2003, hereafter BC03). 

The basic data reduction and extraction of the kinemat- 
ics are described in Paper I. The only change introduced at 
this level is a new criterion for the binning process, following 
the precepts of Cardiel et al. (1998). The aim of this proce- 
dure is to estimate the minimum signal-to-noise ratio (S /N) 
required for the new bins in order to obtain reasonably small 
uncertainties in H/3 line strength measurements (5(11/3) be- 
tween 0.04 and 0.3) and so, in the relative ages. An estimate 
of the minimum 5(11/3) was made for each galaxy by us- 
ing different H/3-versus-metallic-index diagrams and BC03 
model grids; depending on the position of the galaxies' data 
points on the grids, a given uncertainty in the indices will 
translate into a corresponding error in age. In the majority 
of cases, an uncertainty no larger than 4 Gyr was allowed 
for the outermost bins. Then, the following expression of 
Cardiel et al. was used to obtain the S/N required for a 
given uncertainty 5(H/3): 

7.301 - 0.2539 • H/3 



S/N [A] 



5{m 



(1) 



The S/N per A of the new bins are 100, 50 and a minimum 
of ~ 30 (this last value varying somewhat from galaxy to 
galaxy), decreasing as the radius increases. As a result of 
the new binning, the data cover median radii ~ 2 bulge 
effective radius (Re) for this sample of SOs. 



2.1 Transformation to the Stellar Library 
Resolution 

The Lick/IDS indices studied are listed in Table [1] To be 
able to compare the results to simple stellar population 
models, the indices must be measured at the same spec- 
tral resolution as the model's stellar libraries. The width of 
the observed spectral lines are a convolution between the 
instrumental width and the velocity dispersion of the stars 
projected along the line-of-sight. Therefore, to transform the 
galaxies' spectra to the required resolution, they were con- 
volved with a Gaussian of dispersion, a, equal to 



(2) 



where crjib is the resolution from the stellar library, ainst is 
the instrumental widening (~ 30kms~^, see Paper I) and 



Table 1. Lick indices used in this study (Worthey et al. 1994; 
Worthey &; Ottaviani 1997) and resolution at which each in- 
dex was measured when Lick resolution was adopted (Sanchez- 
Blazquez 2004). 



Index 


Lick Resol. 




[kms-l] 


H/3 


225 


Fe5015 


200 


Mgi 


200 


Mg2 


200 


Mg6 


200 


Fe5270 


200 


Fe5335 


200 


Fe5406 


200 


Fe5709 


200 


Fe5782 


200 



(jgai is the velocity dispersion of the galaxy at the corre- 
sponding radius, already calculated during the extraction of 
the kinematics. 

Once the spectra were broadened, line-strength indices 
were measured in the central regions of the galaxies (within 
Ra/8 of the bulge) using the software INDEXF, developed by 
one of us (NCq). Following the Lick/IDS index definition, 
the program calculates a pseudocontinuum (a local contin- 
uum level) for each spectral feature defined by the means 
within two pseudocontinuum bands located one at each side 
of the spectral feature (see Worthey et al. 1994 and Worthey 
& Ottaviani 1997). Then, the line index is measured with 
respect of the pseudocontinuum by integrating within the 
feature central band 



Index — 



1-^ldA, 



(3) 



where Ad and Ac2 are the initial and final wavelengths of 
the central band, Sa is the fiux of the spectrum at a certain 
A and Ca is the corresponding pseudocontinuum fiux. This 
program also estimates the uncertainties resulting from the 
propagation of random errors and from the effect of uncer- 
tainties on radial velocity by performing Monte-Carlo simu- 
lations. 

The choice of BC03 models instead of other alternatives, 
like the models of Thomas, Maraston & Bender (2003, here- 
after TMB03), was made based on the higher resolution of 
the stellar libraries of the former (3 A, Le Borgne et al. 2003), 
which permitted the study of the spectral features of these 
galaxies in greater detail. For the study of the main stel- 
lar population properties, however, a BC03 model with Lick 
resolution had to be used (resolution from FWHM ~ llA 
around 4000 A to ~ 8.5 A around 5000 A) given the higher 
reliability of these models when individual ages and metal- 
licities are estimatecjf]. We emphasise that at both. Lick and 



^ http:/ /www. ucm.es/info/Astrof/software/indexf/indexf.html 
^ The Lick resolution models make use of the so-called fitting 
functions (Worthey et al. 1994; Worthey & Ottaviani 1997) while 
the 3 A models don't. As a result, the model grids of the later 
sometimes show irregular (unphysical) patters which translate in 
very uncertain individual ages and [Fe/H]. 



so galaxies in Fornax: masses and central stellar populations 3 



3 A resolutions, we are using directly the predictions pro- 
vided from BC03. The normal distribution of BC03 mod- 
els includes predictions for both spectral resolutions, so we 
refer the reader to BC03 paper for further details concern- 
ing these models. In synthesis, all indices were measured at 
both resolutions by applying the previously described proce- 
dures. Further comparisons between the two index sets are 
presented in the following sections. 

For NGC 1316, NGC 1380 and NGC 1381, the velocity 
dispersions within their central regions are higher than the 
3 A stellar library's resolution. Therefore, the procedure de- 
scribed above cannot be applied and the indices must be 
corrected after they have been measured. In these cases, we 
proceed as follows: for each galaxy, the best stellar-template 
combination for the affected bins were used as a model of the 
non-convolved galaxy spectra. Using equation [5] these spec- 
tra were widened to the stellar library resolution and then, in 
steps of 20kms~^, convolved with Gaussians of dispersions 
Cgai between and 400kms~^. The next step consisted of 
measuring all line indices in each spectrum and calculate a 
correction factor for each index, C((T)indcx, of the form 

Index — 

lndex(0)/lndex(cr). (4) 

Index(O) is the index measured from the best stellar- 
template combination degraded at the stellar library res- 
olution, and Index((T) is the index measured from the same 
spectrum but widened by crkms"^. Therefore, the corrected 
index is given by 

Index(O) = C(cr) Index ■ Index(cr). (5) 

Notice, however, that for the molecular indices Mgi and 
Mg2, measured in magnitudes, the correction applied was 
of the form 

C(cr) = Index(O) - Index(cr). (6) 

The correction factors for the individual bins turned out to 
be similar within each galaxy. Therefore, it was decided to 
apply a single set of correction factors for each SO by taking 
the average of the individual bins' corrections and fitting the 
results with polynomials of order 3. For NGC 1316 and 1380, 
analogous corrections were applied when Lick resolution was 
used. 

The corrections were applied along the radius, R, to all 
spectral bins for which the velocity dispersion was higher 
than the stellar library's resolution (cr|ai(^) + ''"hist ^ ""ub)- 
For NGC 1380 and NGC 1381 the different corrections 
ranged between 2-10% of the value of the measured in- 
dices, while for NGC 1316 they were slightly larger, typ- 
ically between 10-20% of the original measurements. For 
each galaxy, an estimate of the uncertainty in these correc- 
tions was attempted by measuring the deviations of all the 
coefficients C((j)indox from each bin at a velocity dispersion 
of 400kms~^. The resulting errors translate into uncertain- 
ties of the order of 0.001 A for the different indices. Such 
errors are negligible compared to other sources of error from 
the data reduction/kinematics extraction processes, and so 
were neglected in subsequent analysis. 



2.2 Emission and Lick Spectrophotometric 
Corrections 

For a long time early-type galaxies were considered as 
gas/dust free objects. Ifowever, subsequent work on large 
samples of ellipticals has revealed that about 50% of these 
objects show weak nebular emission lines in their optical 
spectra (Caldwell et al. 1984; Phillips et al. 1986; Goud- 
frooij et al. 1994). Measurements of different nebular lines 
indicates the presence of lO'^-lO^ M© of ionised gas in the 
central regions of these galaxies. Despite the small amount 
of gas, some absorption line features can be affected by these 
emissions. 

Weak traces of nebular emission ([OiiiIasoot rest frame) 
were found in the spectra, which affects the measurement 
of the indices Fe5015 (within its central band) and Mgfe 
(within one of the pseudocontinuum bands). To correct these 
indices of emission, the best stellar-template combination 
was compared to the galaxy spectrum and the spectral fea- 
ture of [Oiii]a5007 was replaced by the corresponding section 
of the best stellar-template combination before the indices 
were calculated. No other signatures of emission, such as 
[Oni]A4959 or [NI]a5200, were found which could affect the 
measurements of other metallic indices. 

It was also found that the index H/3 has traces of con- 
tamination from emission. The correction of this index is 
of particular importance in this study given that it is the 
main age indicator. The presence of emission in H/3 de- 
creases the magnitude of the measured index and the in- 
ferred age becomes older (H/3 decreases with the age of the 
stellar population). Gonzalez (1993), studying a sample of 
bright elliptical galaxies, found an empirical correlation be- 
tween the EW[Oni]A5007 and the EW of H/3 in emission. 
For his brightest galaxies he found that ^^fl^^T^""^ — ~ 0.7. 

=■ ° lUnilA5007 

However, Trager et al. (2000) found in a sample of 27 ellip- 
tical galaxies that this relation has variations between 0.33 
and 1.25, and that the mean value is 0.6 instead of 0.7. Ac- 
cordingly, it was decided to apply the following (additive) 
emission correction for the index H/3: 

A(H/3) = 0.6 ■ EW[Oiii]a5007, (7) 

where EW[Oiii]a5007 was measured from the residual spec- 
tra obtained by subtracting the best combinations of stellar 
templates from the galaxy spectra. The variations found by 
Trager et al. (2000) in the ratio H/?/[Oiii]a5007 are trans- 
lated into uncertainties in the determination of the age of 
the order of 3%, negligible compared to other sources of er- 
ror for this index. For the SO galaxies analyzed here, the 
corrections applied to H/3 were not usually larger than 10% 
of the original measurements. 

When models based on Lick/IDS libraries (e.g. some of 
the models published by BC03, TMB03) are used, a spec- 
trophotometric correction is usually applied to the measured 
indices. This is because the Lick libraries are not calibrated 
in fiux, but flux-normalised by using a tungsten lamp. The 
indices are usually corrected by additive factors calculated 
from observations of stars from the Lick library with the 
same instrumental configuration and photometric conditions 
as the data. We lack such data in the current study; how- 
ever, there are other ways around this issue. We attempt to 
calibrate our Lick resolution indices by using data from the 
literature. In his sample of early-type galaxies, Kuntschner 
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(2000) presents lower resolution, fully calibrated central line 
indices for our entire sample of SOs, being an ideal dataset for 
our purpose. All our indices but Fe5782 were also measured 
in his work. We rebinned our data by using a central aper- 
ture of 3.85" as in Kuntschner (2000) in order to perform the 
comparison. By using the individual measurements for each 
galaxy, mean corrections (offsets) were estimated for each 
index. The uncertainties of these corrections, however, were 
quite large. After studying the significance of the offsets by 
using a Student-t correlation test, we conclude that for all 
the indices but Mgi and Mg2 the offsets were not signif- 
icant at 95% confidence level. These two molecular indices 
are the most likely to be affected by spectrophotometric cal- 
ibration problems given their highly separated pseudocon- 
tinuum bands. In any case, the corrections estimated for the 
remining indices were usually fare below 10% the value of 
the original measurements. Therefore, we only applied the 
corrections to the Mgi and Mg2 indices measured at Lick 
resolution; 



5(Mgi)Lick = +0.0166 ± 0.0056 mag 
5(Mg2)Lick = +0.0216 ± 0.0072 mag. 



(8) 
(9) 



We stress that these corrections were applied only to Lick 
resolution measurements. The stellar libraries on which 3 A 
resolution models are based have being properly calibrated, 
so no corrections are required for 3 A resolution data. 

In this work, two combined indices were also used, (Fe) 
(Gorgas, Efstathiou & Aragon-Salamanca 1990) defined as 

,^ , Fe5270 + Fe5335 

(Fe) = , (10) 

and [MgFe]' (Gonzalez 1993; Thomas, Maraston & Bender 
2003) defined as 

[MgFe]' = ^Mgb • (0.72 x Fe5270 + 0.28 x Fe5335). (11) 

[MgFe]' will be particularly important for the analysis 
of stellar population parameters, like age and metallicity, 
given its almost null dependency on a-element abundance 
(TMB03). All measured and corrected central indices for the 
Fornax sample are presented in Tables lAll fat 3 A resolution) 
and IA2I (at Lick resolution) in the Appendix. 

Finally, note that in different parts of this work the 
atomic indices are sometimes expressed in magnitudes (fol- 
lowing CoUess et al. 1999) according to the expression 



Index* 



-2.5 ■ 



log(l 



Index \ 



(12) 



were Index* represents the index in magnitudes. Index is the 
index expressed in A, and AA is the central bandwidth in A, 
listed in Table 1 of Worthey et al. (1994). In the particular 
case of (Fe), the index in magnitudes was calculated using 

Fe5270 Fe5335 



(Fe>* = -2.5-log- 1-— — , (13) 

Z \ AAfo5270 AAfc5335 

while [MgFe]' was always expressed in A. 



3 COMPARISON BETWEEN INDICES AT 
DIFFERENT RESOLUTIONS 

In this section, comparisons between line indices measured 
at different resolutions are made in order to understand in 



a qualitative way the possible effects on the ages and metal- 
licities when inferred from degraded spectra. In Figure [TJ 
measurements of all central indices at 3 A and Lick resolu- 
tion are compared for the whole galaxy sample. 

All forthcoming comparisons are made in relation to 
the 3 A indices. The H/3 index suffers little degradation, 
where only the galaxy NGC 1375, with the strongest H/3- 
absorption, may suffer from a small relative decrement of 
its luminosity-weighted age estimate. The Mg indices suf- 
fer a dissimilar behaviour. While Mg6 presents a very little 
variation between the two resolutions, the molecular indices 
Mgi and Mg2 show a systematic offset with respect to the 
1:1 lines. Their larger values at Lick resolution are just a 
consequence of applying the spectrophotometric correction 
in Section[2]2] In general, Fe indices present systematic shifts 
with respect to the 1:1 lines with slightly larger differences 
for galaxies with strong spectral features. The dependence 
on index intensity may introduce some differences in the rel- 
ative metallicities derived using Fe features by diminishing 
this parameter in galaxies with strong Fe lines. However, 
the relative differences in shifts among the galaxies are not 
strong enough to make this an important effect. The relative 
index variations between galaxies can be more easily appre- 
ciated in Figure (2] where the indices at 3 A resolution are 
plotted versus the fractional change when they are measured 
at Lick resolution. 

The Fe features vary, on average, by —15%, while the 
H/9 and the Mg& index variations are usually less than —10%. 
The corresponding panels for Mgi and Mg2 ilustrate not 
only the systematic offsets of Figure [T] but the magnitude 
of the spectrophotometric corrections applied to these two 
indices. The variations are particularlly important for Mgi, 
where a median percentual change of +20% is found. This is 
the main reason why we are going to restrict the use of Mgi 
and Mg2 to the qualitative tests presented in the following 
sections, without attempting to estimated individual ages 
and [Fe/H] with them. In this way, Mgb becomes our more 
reliable magnesium tracer. Concerning the estimate of the 
[Mg/Fe] (or [a/Fe]) relative abundance tracers, the fact that 
the Mg& index presents little variation in comparison to the 
Fe indices will certainly introduce changes for each individ- 
ual galaxy if the abundance tracers are measured as the dif- 
ferences between metallicities (e.g. [Fe/II]Mgb — [Fe/H](Fc))- 
However, the similar variations in the Fe indices experienced 
by all the galaxies will result in the relative [Mg/Fe] abun- 
dances between them being mostly unaffected. In any case, 
the index ratio Mg&/(Fe) provides an alternative model- 
independent tracer of the a-element abundances previously 
applied by other authors (see TMB03). 

In summary, in this section comparisons between line 
indices measured at 3 A and Lick resolution have been made 
in order to find possible systematic differences which could 
affect our conclusions when one set or the other is used. 
Negligible variations are found in the H/3 and Mg6 indices. 
A somewhat larger change is detected in the Fe features, but 
we conclude that such effects will not introduce important 
relative changes in the ages, metallicities and abundance 
tracers between the galaxies. 
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Figure 1. Comparison of the central line indices measured at 3 A and Lick resolutions. 



4 THE MODELS: CONSISTENCY WITH THE 
DATA 

Before sliowing tlie main results, in this section some con- 
sistency clieclis between the BC03 models and the data are 
presented. These tests give an idea of how well the models 
describe the data and, in consequence, how the model results 
should be interpreted. For a stellar population study, it is 
important to check if there are significant differences in the 
relative spectrophotometric calibrations of the data and the 
stellar libraries used by the models. If differences are found, 
the results must be interpreted accordingly. Although ages 
and metallicities were estimated using Lick resolution data 
and models, some qualitative comparisons were also made 
with the higher resolution data. Therefore, the consistency 
tests were performed for both the 3 A and Lick resolution 
BC03 models and data sets. 



4.1 3 A resolution models and data 

When using the 3 A resolution STELIB library, it is expected 
that little difference from the data will be found because 
both sets of spectra were flux calibrated using spectropho- 
tometric standard stars. 

This assumption was tested by creating index-index 
plots which are almost degenerate in age and metallicity. 
In this way, the model grids describe narrow bands in the 



index-index space, regions which should trace the galaxies' 
data if the models describe accurately the properties of these 
objects and if there are no problems with the relative flux 
calibrations. 

In Figure [S] different Mg indices are plotted against 
each other. The folded grids represent BC03 3 A resolution 
models, and the line indices measured within J?o/8 are plot- 
ted for the complete SO sample. The match between data 
and models is very good, with a hint of a small systematic 
deviation for galaxies with intense Mg indices (NGC 1316, 
1380 and 1381). Even less pronounced, a small offset to- 
wards lower Mgl values is apparent from the Mgi versus 
Mgb panel. A similar plot is presented in Figure |4] for six Fe 
indices. In these cases, the scatter from the models is larger 
than for the Mg indices, but the models seem to trace the 
mean trends in the data. For the three plots with Fe5015, 
however, there seems to be some systematic deviation be- 
tween the data and model grids. Because the central band- 
pass of Fe5015 was corrected from [Oiii]a5007 emission, this 
particular index must be used with caution. 

Thus, as predicted, no important mismatches were 
found between indices and model predictions at 3 A resolu- 
tion, suggesting that the relative flux calibrations are consis- 
tent with each other. By assuming solar abundances, BC03 
models are able to consistently reproduce the observations 
for the majority of these SO galaxies. Therefore, in the fol- 



6 A.G. Bedregal, A. Arag on- Salamanca, M.R. Merrifield and N. Cardiel 



o 

a 

A 




2 3 4 

H|S (31) [A] 




0.05 0,1 0.15 
Mgl (3A) [mag] 




0,1 0.15 O.S 0.H5 o.a 
MgS (3i) [mag] 



tj 

3 

A 

a 



12 


L 


1,11 


1 1 , 1 w 1 , 1 1 1 1 ,J 


14 






• J 


le 
le 






' •'• ^ 


so 








22 


r, , 


• 

1 1 1 1 


, , 1 , , , , 1 , , 




3 


4 


5 6 7 






Pesois (aA) [A] 


10 
13 


il ' 


' 1 


•k' ' ' L 

« 

• 


14 










-I 1 1 




• 

. < . 1 . . < . 1 . . . . 1- 




1.5 


3 


3,5 3 3,E 






re6370 (aA) [A] 


14 


.11,1 






le 






• 


IB 








30 








S3 


-nil 


llffi 


- 




1,5 


s 


E,5 3 a, 5 



re53a5 (aA) [A] 




0,6 o,a 1 1,3 
Fesvoa (3A) [A] 




0.4 O.S 0.8 1 
Fe57B3 (aA) [A] 




1,5 a 3,5 a a, 5 

<¥e> (3A) [A] 



1- 1 1 , ,^ , , 1 1 , , , , , 1 , , ,- 
: • • 

- • • - 


-13 
-14 


: 1 • 




1 

H 


• - 


-S 
-10 


. 1 1 1 1 


* 

• • 

* 


* .• 


-le 

-18 










-IS 


• 




1 • ■ 




", 1 , , 





m 




-14 


-, , 1 , 


• 1 , ,f 


3 3 4 5 




1 


1.5 


a 


3.5 


3 


3 4 



Mgb {3A) [1] 



re6406 (aA) [A] 



[MgFe]' (3A) [A] 



Figure 2. Central line indices at 3 A resolution versus fractional change (%) when measured at the Lick resolution. 



Fe5015, 





o 

^ 1.5 



o 



(3l) 



[A] 





1 1 1 1 laJ 1 I// 







2 2.5 3 3.5 0.6 0.8 1 

<Fe>,3i) [A] Fe5709,3i. [A] 




Figure 3. Consistency test using Mg indices and BC03 models at 

3 A resolution. Folded grids arc BC03 models and datapoints cor- 
respond to central indices (within Rg/S) of SO galaxies in Fornax 
Cluster. 



Figure 4. Consistency test using Fe indices and BC03 models at 
3 A resolution. Folded grids are BC03 models and datapoints cor- 
respond to central indices (within Re/S) of SO galaxies in Fornax 
Cluster. 
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Figure 5. Consistency test using Mg indices and BC03 models at 
the Lick resolution. Folded grids are BC03 models and datapoints 
correspond to central indices (within Rc/8) of SO galaxies in the 
Fornax Cluster. 



Figure[6]shows Fe indices in good agreement with model 
predictions. The scatter found in the previous section and a 
hint of a systematic deviation in panels involving Fe5015 are 
also present here, which reinforces the idea that the Fe5015 
index must be used with caution. After the spectrophoto- 
metric correction applied in section 12.21 a very good match 
is also found between Mg indices and model grids. In his 
study of early type galaxies in Fornax, Kuntschner (2000, 
hereafter KOO) found in similar tests discrepancies between 
the Mg indices and other Lick-based models (Worthey et al. 
1994 and Vazdekis et al. 1996). We do not find such differ- 
ences here. We cannot totally discard, however, the existence 
of such systematic discrepancies found by KOO given (i) the 
large uncertainties of our Lick spectrophotometric correc- 
tions (derived from a small sample) and (ii) the relatively 
large errors of our Mgi indices. 

We cannot discard other sources for the individual off- 
sets between models and data (at both. Lick and 3 A res- 
olutions). For example, the index Mgi is highly dependent 
on the element carbon, which could explain some deviations 
observed in Figures [3] and |S] if some peculiarities in carbon 
are present. 

In any case, given the good agreement between models 
and data, and having in mind the discussion in Section [2.11 
(see footnote ^), we can confidently use these models to cal- 
culate ages and metallicities for individual galaxies. 



Fe5015(,,,,, [A] 








r 




Figure 6. Consistency tost using Fe indices and BC03 models at 
the Lick resolution. Folded grids are BC03 models and datapoints 
correspond to central indices (within Re/ 8) of SO galaxies in the 
Fornax Cluster. 



lowing sections these models can be confidently used to make 
qualitative comparisons with this data set. 



4.2 Lick resolution models and data 

The test described in the previous section were also applied 
to the Lick resolution data and models and the results are 
presented in Figures [5] and [6] 



5 RESULTS AND DISCUSSION 

In this section, the main results of the central (within Ra/8) 
stellar populations study are presented. As mentioned ear- 
lier, KOO studied the central stellar populations of all ellipti- 
cals and SOs in the Fornax Cluster using lower quality optical 
spectra than the ones presented in this work. Therefore, his 
study will be a constant reference in the forthcoming pages. 

From this point on, the BC03 models at 3 A resolu- 
tion were always preferred for qualitative comparisons be- 
tween the data and model grids. Also, for the correlations 
between line indices versus different dynamical mass tracers 
(e.g. Index*-log(cro) relations) the 3 A resolution data were 
used in order to study the spectral features in grater detail. 
Only for the explicit calculation of ages, metallicities and 
abundances tracers (e.g. [Fe/H]Mgb — [Fe/H] (pc) , Mg6/(Fe)) 
were the models and data at the Lick resolution applied. Un- 
less explicitly mentioned, the 3 A resolution data and models 
should be assumed. 

Following previous works in this field (e.g. KOO; Mehlert 
et al. 2003; Sanchez-Blazquez et al. 2006), we focus our study 
of central stellar populations on the relations between dif- 
ferent line indices and kinematical parameters of the host 
galaxies. The relations between 10 line indices and the cen- 
tral velocity dispersion, ao, are studied in order to unmask 
the driving variables behind the observed slopes and scatter. 
In Figure [T] the relations between line indices and log((Jo) 
are presented for the sample of 9 SOs. 

Rather than forming continuous trends. Figure [7] shows 
the three bright galaxies (NGC 1316, 1380 and 1381) and 
the rest of faint SOs in two separate clumps, each one in 
opposite extremes of the ao range. However, because the 
sample was selected only according to morphological classi- 
fication at optical wavelengths (see KOO), there is no obvi- 
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Figure 7. Central line indices (in magnitudes) versus log(cro). Continuous lines are the best fits to the datapoints; red dashed lines 
represent the best fits to normal ellipticals in Fornax from Kuntschner (2000). 



ous selection bias against galaxies of intermediate ao values. 
Therefore, this sample will be considered to describe the two 
extremes of continously-populated Index* -log (cro) relations. 

The plots in Figure [7] reconfirm some of the results al- 
ready indicated by KOO. First, faint SOs have a large scatter 
with respect to the main relations, while the bright galaxies 
NGC 1380 and NGC 1381 closely follow the relations for 
ellipticals estimated by KOO (red dashed lines in Figure [7]). 

Second, the merger remnant NGC 1316 (triangle in Fig- 
ure[7)l departs significantly from the Mgi- and Mg2-log((To) 
relations (at more than 4-sigma level in both indices) to- 
wards low Mg absorption at a given ctq, while closely fol- 
lowing the remaining Fe-log((Jo) relations. These facts, and 
the comparatively high H/3 absorption in the H/3-log((To) 
diagram have been interpreted by KOO as the signature of 
the presence of a younger stellar population in this galaxy, 
probably formed during the merger. Such a young popu- 
lation could reduce the intensity of some metallic features 
by increasing the overall continuum level. This effect would 



be particularly important for strong metallic features, like 
the Mg triplet (indices Mgb and Mg2) and for the molec- 
ular indices Mgi and Mg2 whose pseudocontinuum level is 
strongly dependent on large-scale changes in the continuum. 
The weaker Fe features in this spectral range could be less 
influenced by this effect, explaining why NGC 1316 appears 
in good agreement with the mean Fe relations. 

Finally, the slopes of the Mg-log((To) trends are steeper 
than the ones for the Fe lines (see fit equations in Table [5}, 
as KOO also found for his fits to normal ellipticals. 



5.1 The Slopes of the Index*— log((7o) relations 

To study which parameters are driving the observed slopes 
in the Index*-log(cro) relations, the models of BC03 were 
used to parametrise the individual indices as a function of 
metallicity ([Fe/H]) and log(age) such that 

IndexBC03 = A + B- [Fe/H] + C ■ log(age[yr]), (14) 
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Table 2. Parameters of the linear fits Index* = a + h ■ log((To) of 
the SO galaxies in Fornax, a^^^ is the standard deviation about 
the fit. 



Mgb„,,y [1] 



Index 


h±5{b) 


a it &{a) 


0"std 




H/3* 


-0.0633 ± 0.0031 


0.2114 ± 0.0066 


0.0360 




Fe5015* 


0.0107 ± 0.0021 


0.0625 ± 0.0045 


0.0219 




Mgi* 


0.1290 ± 0.0016 


-0.1618 ± 0.0034 


0.0229 




Mg* 


0.1683 ± 0.0018 


-0.1111 ± 0.0038 


0.0462 




Mgb* 


0.0932 ± 0.0027 


-0.0486 ± 0.0056 


0.0258 




Fe5270* 


0.0185 ± 0.0022 


0.0509 ± 0.0045 


0.0161 




Fe5335* 


0.0303 ± 0.0025 


0.0262 ± 0.0052 


0.0135 




Fe5406* 


0.0293 ± 0.0027 


0.0214 ± 0.0056 


0.0112 




Fe5709* 


0.0077 ± 0.0025 


0.0349 ± 0.0053 


0.0125 




Fe5782* 


0.0284 ± 0.0032 


-0.0113 ± 0.0068 


0.0071 


S 


(Fe>* 


0.0538 ± 0.0037 


0.8272 ± 0.0077 


0.0310 


[MgFe]' 


1.5212 ± 0.0801 


0.4087 ± 0.1661 


0.5311 





where IndexBQQ3 are index values from BC03 models and 
9(Index^co3) 



B = 



C = 



a[Fe/H] 

^(IndexBcoa) 
a(log(age)) 



(15) 
(16) 



Because the galaxies in this sample cover a wide range 
in ages and [Fe/H], we use models with ages from 1.0 to 
12.6 Gyr and metallicities from —0.64 to 0.55 dex. The re- 
sulting parametrisations are presented in Table [S] The last 
two columns give the required dependence of log((To) on 
[Fe/H] and age respectively, if that is the only driving vari- 
able for the observed slopes in the Index*-log((To) relations. 
This dependency was estimated for each index by combining 
the parametrisations from the models (egs. I15l and ll6p with 
the slopes measured in Figure [7] (6 in Tabled. If [Fe/H] 
is considered as the driver of the observed Index*-log((To) 
trends 



b_ _ a(Index'') a[Fe/H] 
B a(log(cro)) 9(Index*co3) 



a[Fe/H] 



(17) 



(18) 



9(log(cro)) ' 

Alternatively, if age is taken as driving variable 

b 9(Index*) cl(log(age)) 9(log(age)) 
C a(log(cro)) 9(Index*co3) a(log(cro)) ' 

Interestingly, these parametrisations show that the Mg- 
log((To) relations depend as strongly on age as on metallic- 
ity. On the other hand, the slopes of the Fe indices versus 
log((To) relations can be reproduced with weaker dependen- 
cies on metallicity than on age, but the differences between 
both dependencies are not very large. Naively, a much lower 
dependence on metallicity relative to the one on age would 
be expected for these metallicity tracers. Only the behaviour 
of the H/3*-log((To) relation matches the expectations, pre- 
senting a slope much easier to reproduce by a small age 
dependence than by a huge metallicity one. 

It is also interesting to note that not all the [Fe/H]- 
sensitive indices have the same dependence on [Fe/H]; the 
three Mg indices (Mgi, Mg2 and Mgb) required a stronger 
[Fe/H] dependence than the Fe indices in order to explain 
their slopes in the Index* -log (ao) diagrams. This would not 
be expected if metallicity is the only parameter varying with 
log(cro). 




0.6 8 1 1.2 1.4 1 6 



Fe5709,„„„, [A] 



Figure 8. Central metallic indices versus H/3 for SOs in Fornax 
Cluster. The grids represent BC03 models labeled in terms of age 
and [Fe/H]. 



Therefore, these results suggest that age could be an 
important driver of the observed correlations while the role 
of metallicity is more uncertain. Also a third variable, such 
as the relative Mg abundances, may affect the relations by 
producing relative differences between indices, each one with 
different sensitivities to changes in the chemical composition 
of these systems (Vazdekis et al. 2001; Sanchez-Blazquez et 
al. 2006). 



5.1.1 The [Fe/H]- and log(age)-log(ao) relations 

A closer comparison between the data and the BC03 models 
could shed more light on the origin of the observed cor- 
relations. Individual ages and [Fe/H] were calculated by 
plotting different [Fe/H]-sensitive indices versus H/3; data 
and models at the Lick resolution were used in this step. 
The details of the procedure applied here are described in 
Cardiel et al. (2003). Briefly, the ages and metallicities were 
logarithmically-interpolated from the model grid values by 
fitting the four closest cells to a datapoint with quadratic 
polynomials. The procedure also considers the effect of co- 
variance when estimating the uncertainties for both param- 
eters. In this way the age-[Fe/H] degeneracy is broken and 
luminosity-weighted estimates of both parameters were ob- 
tained. As an example, in Figure [S] four of these plots are 
presented for the central index values, also including the 
BC03 model grids labeled in terms of age and [Fe/H]. Cen- 
tral ages and metallicities for different indices are presented 
in Table IX3l for all the Fornax sample. 

In Figure [9l log(cro) is plotted against age and [Fe/H] 
using different [Fe/H] sensitive indices. The continuous lines 
represent the best linear fits whose slopes are presented 
in Table |4] It is surprising that no obvious trends are ob- 
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Table 3. Parametrisation of the line indices using BC03 models: Indexg^^^g = A + B ■ [Fe/H] + C ■ log(age[yr]). Last two columns: 
dependence of log(a'o) on metallicity (age) if this is the only driving variable for the observed slopes in the Index*-log(a"o) relations. 



Index 


A 


B±5{B) 


C ± 5{C) 


a[FG/H] 

8(log{CTo)) 


e{log(agG)) 

a{iog(o-o)) 


H/3* 


0.7778 


-0.0157 ± 0.0014 


-0.0712 ± 0.0020 


4.029 


0.888 


Fe5015* 


-0.0854 


0.0323 ± 0.0007 


0.0177 ± 0.0011 


0.330 


0.601 


Mg* 


-0.4602 


0.0502 ± 0.0023 


0.0547 ± 0.0032 


2.568 


2.357 


Mg* 


-0.8310 


0.1075 ± 0.0031 


0.1053 ± 0.0044 


1.564 


1.597 


Mg6* 


-0.4431 


0.0521 ± 0.0014 


0.0572 ± 0.0019 


1.790 


1.629 


Fe5270* 


-0.1509 


0.0366 ± 0.0005 


0.0244 ± 0.0007 


0.505 


0.757 


Fe5335* 


-0.1665 


0.0388 ± 0.0006 


0.0256 ± 0.0008 


0.779 


1.180 


Fe5406* 


-0.2017 


0.0419 ± 0.0006 


0.0284 ± 0.0008 


0.699 


1.029 


Fe5709* 


-0.0511 


0.0199 ± 0.0007 


0.0098 ± 0.0009 


0.386 


0.781 


Fe5782* 


-0.0856 


0.0293 ± 0.0006 


0.0137 ± 0.0009 


0.968 


2.074 



Table 4. Slopes from linear parametrisation of log(age) or [Fe/H] 
versus log(o'o) for SO galaxies in Fornax (Figure|9]l using different 
line indices as metallicity tracers. 



Parameter 


Slope 


Parameter 


Slope 


log{age)[MgFo]' 


0.946±0.079 


[Fc/H][MgFcl' 


0.208±0.090 


log{age)(po) 


0.994±0.075 


[Fe/H]<p,) 


-0.189±0.057 


log(age)Mgb 


0.803±0.080 


[Fe/H]Mgb 


0.806±0.106 


log(age)Fe5709 


1.033±0.074 


[Fe/H]Fo5709 


-0.189±0.113 



served in any of the [Fe/H] -log (cro) plots, considering that 
the Index-log ((Jo) relations are usually attributed, at least 
for ellipticals, to a chemical composition change with galaxy 
mass. The slope for the [Fe/H] -log (ctq) relation using Mg& 
seems to be an artifact of the linear fit instead of pointing to 
a real correlation. Comparing the slopes for metallicity pre- 
sented in Tabled with the required metallicity dependences 
for each Index*-log(cro) relation in Table |3] it is clear that 
[Fe/H] cannot be the main driver of any of the Index* corre- 
lations with log((To). The plots against age present stronger 
trends in the same direction as the Index*-log((To) diagrams: 
the three high- dispersion galaxies seem to be, on average, 
older than the rest of the sample. This result holds for all 
the age estimations, including the one using [MgFe] ' which is 
insensitive to a-element relative abundance variations. Com- 
paring the slopes of these relations with the required age de- 
pendences for the observed Index* -log (ctq) trends, we found 
that the H/3* and almost all the Fe trends can be explained 
by an age effect. However, differences in relative age do not 
seem to explain the steeper slopes found for the Mg indices 
in comparison to Fe, reinforcing the idea that a-element rel- 
ative abundances are playing an important role too. 



5.1.2 The [Mg/FeJ-log(ao) relations 

At the beginning of this section, the possibility of a Mg (or 
Q-element) abundance dependence of the Index*-log((To) re- 
lations was suggested, and this possibility was re-enforced by 
the results presented above. This alternative was further ex- 
plored by plotting Mg-sensitive versus Fe-sensitive indices 
and comparing the data to the solar-abundance predictions 
of the models (see Worthey, Faber & Gonzalez 1992). In 
Figure IIOI twelve panels plotting central Mg versus Fe in- 
dices are presented, together with BC03 model grids. What 
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Figure 9. Central Ages and metallicities versus log((To) for For- 
nax SOs. Lines represent the best linear fits to the data. 



is clear from almost all these panels is that the three highest- 
dispersion galaxies in the sample (open symbols in Fig- 
ure [TUl) are not accurately described by models that assume 
solar abundance ratios, presenting stronger Mg indices for 
a given Fe index than the model predictions. Among these 
three galaxies, the effect seems to be stronger for NGC 1380 
and 1381, while NGC 1316 (the merger remnant) is usually 
closer to the solar-model grids. In similar plots, KOO found 
that only the nucleus of NGC 1380 presented this behaviour, 
suggesting that this galaxy is Mg overabundant compared 
to Fe. 

In order to trace the relative Mg overabundances we 
compare the relative metallicities derived from Fe- and Mg- 
sensitive indices and the same age indicator (H/3 in this 
case). At this point, it is worthy to stress that we are not at- 
tempting to derive absolute Mg (or a-elements) abundances. 
By comparing different metallicity estimates (and by using 
Mg6/{Fe), see below) we aim to trace the relative Mg abun- 
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Figure 10. Central Mg versus Fe indices for Fornax SOs. Grids correspond to BC03 models at 3 A resolution. 



dances. AH our results and conclusions are independent of 
absolute abundance ratios. 

In the top panel of Figure fTTl \og{ao) is plotted against 
the difference of the metallicities calculated by using Mg6 
and (Fe) as metallicity indicators. Data and models at the 
Lick resolution were used to estimate these metallicities. 
It is apparent from this plot that bright SOs are, on av- 
erage, Mg overabundant relative to fainter galaxies; how- 
ever the uncertainties coming from the models make this 
difference only marginally significant. The bottom panel of 
Figure [TT1 shows the index quotient Mg6/{Fe) versus the ve- 
locity dispersion. The ratio of these two indices has been 
used by TMB03 to calibrate the [a/Fe] abundance ratio of 
their models, so it traces the relative behaviour of the Mg 
overabundance and, equally important, it provides a model- 
independent test. With much reduced errors, it is clear that 
bright SOs have larger Mg&/(Fe) ratios than fainter objects. 
This model-independent test can explain why the slopes of 
the Index*-log((To) correlations are steeper for the Mg in- 
dices than for the Fe ones, confirming the previous results 



of this section. The overall results also imply the existence of 
a correlation between age and Mg&/(Fe) (or [a/Fe]; see Fig- 
ure [TT] bottom panel in first column), in agreement to the 
previous findings of Fisher, Franx & lUingworth (1995) and 
Thomas, Maraston & Bender (2002) in elliptical galaxies. 

Different authors have tried to explain the observed 
[Mg/Fe] overabundances in the context of galaxy formation 
and evolution. All these scenarios are based on nucleosyn- 
thesis models which predict that Type II supernovae mainly 
eject Qf-elements, like Mg, while Type la supernovae are the 
main sources of enrichment of the ISM with Fe. Possible 
explanations include variations in the initial mass function 
(e.g. Schmidt 1963; Worthey et al. 1992; Vazdekis et al. 1999; 
Nagashima et al. 2005), selective loss of metals (Worthey et 
al. 1992) , different time scales of star formation and different 
star formation histories. The last two have been very pop- 
ular in the last few years in order to explain the observed 
overabundances in bright ellipticals (e.g. Mehlert et al. 2003; 
Sanchez-Blazquez et al. 2006). In this scenario, early-type 
galaxies are believed to be remnants of gaseous mergers. 
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Figure 11. Mg/Fe relative abundance tracers versus log((To) for 
Fornax SOs. 

This process produces a starburst where part of the gas is 
consumed and the remainder is ejected via stellar winds (for 
low mass remnants) or heated by different mechanisms (e.g. 
AGN, SNe, thermal conduction) in more massive systems, 
stopping further star formation. According to the hierarchi- 
cal merging paradigm, the starburst is more intense in the 
central regions of the remnant given the stronger potential 
well in the centre and subsequent dissipation of the gas to- 
wards it. The merger process implies a faster build-up of the 
system total mass, having shorter star formation timescales 
than systems grown by secular or more gentle processes. A 
short burst of star formation would not give enough time to 
incorporate the gas enriched with Fe from SNe Type la to 
the formation of further generations of stars, while SNe Type 
II would succeed in incorporation a-elements (like Mg) given 
their much quicker occurrence (in less than 1 Gyr) . In conse- 
quence, a relative [Mg/Fe] overabundance will be observed 
in these systems (understood as an underabundance of Fe) 
with respect to galaxies with younger luminosity-weighted 
ages and more extended star formation histories. 

This scenario is in agreement with the results found here 
for the bright SOs, and is not surprising for NGC 1316, a well 
known merger remnant. However, whether this picture can 
be apphed to the entire NGC 1380 and NGC 1381 galaxies 
or only to their nuclei, is something that cannot be tested 
with central line indices alone. 



5.1.3 The Residuals of the Index* -log (ao) relations 

Previous studies have looked for correlations between the 
Index*-log((To) residuals (AIndex*) and different parame- 
ters that characterise the galaxies' structure or their stellar 
populations. Some of these efforts point out the existence 
of a dependence between AIndex* and the relative ages of 
elliptical galaxies (Schweizer et al. 1990; Gonzalez & Gorgas 
1996; Worthey & CoUobert 2003) by comparing the residu- 
als to the index H/3. Hints of similar trends are seen for the 



Table 5. Parameters of linear fits Index* = a + b ■ log(VMAx) of 
SO galaxies in Fornax. (Tatd is the standard deviation about the 
fit. 



Index 


b±S{b) 


a ± (5(a) 


Cstd 


H/3* 


-0.1099 ± 0.0046 


0.3354 ± 0.0108 


0.0214 


Fe5015* 


0.0197 ± 0.0030 


0.0389 ± 0.0071 


0.0143 


Mg* 


0.2015 ± 0.0026 


-0.3665 ± 0.0062 


0.0233 


Mg* 


0.2593 ± 0.0026 


-0.3633 ± 0.0061 


0.0355 


Mgfe* 


0.1446 ± 0.0038 


-0.1906 ± 0.0089 


0.0194 


Fe5270* 


0.0358 ± 0.0031 


0.0062 ± 0.0073 


0.0070 


Fe5335* 


0.0462 ± 0.0036 


-0.0181 ± 0.0083 


0.0094 


Fe5406* 


0.0475 ± 0.0038 


-0.0279 ± 0.0090 


0.0057 


Fe5709* 


0.0185 ± 0.0037 


0.0080 ± 0.0086 


0.0050 


Fe5782* 


0.0377 ± 0.0047 


-0.0399 ± 0.0109 


0.0098 


(Fe>* 


0.0895 ± 0.0053 


0.7313 ± 0.0122 


0.0166 


[MgFe]' 


2.5019 ± 0.1147 


-2.2356 ± 0.2655 


0.2802 



Mg indices of these SO galaxies but they do not seem to be 
statistically significant, while residuals from the Fe relations 
are clearly uncorrelated with H/3. Potential correlations be- 
tween AIndex* and other parameters were explored, includ- 
ing kinematical (ao; maximum -circular- rotational velocity, 
Vmax) and structural (i?c, Rd, bulge-to-total fraction, Sersic 
index) parameters derived in Papers I and II. No statistically 
significant trends were found for any of these combinations. 

5.2 The Index*— log (Vmax) relations and the 
Dynamical Mass 

Given the strong rotational support of many of these SO 
galaxies, we decided to explore the Index*-log(FMAx) rela- 
tions. Figure [12] shows the different line indices versus max- 
imum (circular) rotational velocity, where clear trends ap- 
pear in all the panels. Linear least-squares fits to the data 
were performed and presented as straight lines in each panel. 
The galaxy with the lowest velocity, ESO 359-G002, was ex- 
cluded from the fitting process because the slope and zero- 
points would be strongly affected by this pressure-supported 
object for which Vmax is highly uncertain. The exclusion of 
the other pressure-supported galaxy, NGC 1316, does not 
significantly affect the fitted slopes and zero-points. 

In Table [5] the linear fit coefficients and standard devia- 
tions, (Tstd, are presented. It is interesting to notice that for 
almost all the indices, aatd is smaller than its counterpart 
from the Index*-log(cro) fits (see Table [2}. If we assume that 
mass is the fundamental physical parameter governing the 
properties of these galaxies, the improvement in the fits may 
be interpreted in the sense that Vmax is a better estimator 
of the total dynamical mass than ao, which is not totally 
surprising after highlighting in Paper I the predominant ro- 
tational nature of these systems. As a test, the mass of these 
galaxies was parametrised as a function of the velocity dis- 
persion and of Vmax in order to compare how central indices 
correlate with both possible descriptions. If the galaxies are 
supported by rotation, we can write 

MxVdiAx-Rd, (19) 

where Af is the total dynamical mass and the disk scale- 
length, i?d, was used to characterise the radius of the galax- 
ies because it reaches the regions where the rotation curves 
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Figure 12. Central line indices (in magnitudes) versus log(VMAx)- Solid lines are the best fits to the datapoints, excluding ESO 359- 
G002. The triangle represents NGC 1316, a merger remnant. For the two galaxies with arrows, the dcprojected azimuthal velocity, V^, 
was used instead of Vmax as a lower limit. 



become flat. If the galaxies are mainly supported by velocity 
dispersion, a similar expression can be obtained, 



(20) 



where Rc is the effective radius of the bulge and {<j'^)r„ is the 
average mean-square dispersion inside R^. Other measures 
of scale such as the half-light radius determined in Paper 2 
produce similar results, but for this analysis R^ is adopted 
because it is well constrained by the GIM2D models (Simard 
et al. 2002). 

In Figures [131 and [T4l central line indices are plotted ver- 
sus Rc ■ (cr^)H,e and R^ ■ Fmaxi respectively. Clearly, both pa- 
rameters, a and Vmax, are tracing, to some extent, the mass 
of these systems. However, an important difference arises be- 
tween pressure- and rotation-supported galaxies (open and 
filled dots, respectively): when the mass is estimated using a, 
it is underestimated in rotationally supported systems with 
respect to pressure-supported ones, as it is evident from the 



offset between open and filled dots in Figure lTSI On the other 
hand, Figure [14] shows that if we estimate the mass using 
Vmax, all the galaxies follow a common sequence. In this 
plot, lower limits for the two pressure-supported SOs (using 
V^) are indicated as open dots with arrows. For these two 
galaxies, masses were also estimated from their measured 
cr's by assuming a simple isothermal sphere in hydrostatic 
equilibrium and presented as open stars in Figure 1141 The 
much tighter trends observed in this plot are consistent with 
our hypothesis that Vmax is a much better tracer of the dy- 
namical mass of these systems than a. 

Figure [151 presents central ages and [Fe/H] from BC03 
models and Mg6/(Fe) versus 3 different dynamical mass 
tracers (Vmax, Ro ■ (o-^)rc and Rd ■ V^max)- No strong cor- 
relations are obvious from the central stellar population 
properties plotted against the dynamical mass. However, for 
■ ^MAx mass tracer, a Spearman rank and a Student-t 
correlation tests (Table|6]) indicate that the relation between 
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Figure 13. Central line indices versus log(i?c ■ o"?,, \)i proportional to the total dynamical mass for pressure-supported galaxies. Open 
symbols represent NGC 1316 and ESO 359-G002, two pressure supported galaxies. 



the [a/Fe] tracer and dynamical mass is significant at a 95- 
99% confidence level. To understand the absence of other 
correlations, we explore the possibility that the uncertain- 
ties in BC03 models are washing out the real trends with age 
and metallicity by considering other SSP models. We have 
used an improved version of Vazdekis et al. (1996) models 
(Vazdekis et al. in preparation; hereafter Vaz96) and re- 
calculate all central ages and metallicities ([M/If] in Vaz96 
nomenclature) for our SO sample 0. In Figure [T^ we present 
the resulting ages and [M/H] versus different mass tracers. 
No strong correlations are observed either with the new 
model predictions, weakening the assumption that model 
artifacts are responsible of the absence of trends between 
the stellar population parameters and the dynamical mass. 
An alternative explanation is that the tight correlations ob- 



^ We have used H/3 and [MgFe] ' indices at Lick resolution as our 
age and metallicity indicators, respectively. 



served between the line indices and galaxy mass are not 
governed by the variation of one of these stellar population 
parameters alone, but by a combination of them in some 
correlated way. 

Indeed, Figure [17] suggests that age and a-element 
overabundance are correlated (97.5% confidence level). 
High Qf-element overabundance implies short star-formation 
timescale. This can be understood by considering a simple 
model for the star-formation history of the Fornax SOs in 
which all these galaxies started forming stars at approxi- 
mately the same time and formed stars for a period At. 
Thereafter, star formation ceased and the galaxies evolved 
passively for a time t until the epoch of observation. The 
age derived from the absorption-line indices is ~ t, while the 
Q-element overabundance is anticorrelated with At. In this 
simple model t is anticorrelated with At, and it follows that 
the measured age must be correlated with the a-element 
overabundance, as observed. If this interpretation is cor- 
rect, the most massive SOs have the shortest star-formation 
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Figure 14. Central line indices versus log(R^ - V^^-yr) , proportional to the total dynamical mass for rotationally-supported galaxies. Open 
symbols with arrows represent NGC 1316 and ESO 359-G002, two pressure supported galaxies. For these two objects the deprojected 
azimuthal velocity, V^, was used as a lower limit instead of Vmax- Open stars correspond to the dynamical mass of these two galaxies 
estimated from a by assuming an isothermal sphere in hydrostatic equilibrium. 



Table 6. Spearman and Student-t tests for central stellar pop- 
ulation parameters versus dynamical mass of SO galaxies in For- 
nax (as estimated from ■ V^^j^, third column in Figure [15}. 
In parenthesis, probability of rejecting the null hypothesis of no 
correlation. 



Dyn. Mass vs. 


Spearman Test Coeff. 


Student-t test Coeff. 


log(age)[MgFcl' 

[Fe/H][MgFo]' 

Mg6/(Fe) 


0.4167 {< 90%) 
0.2833 {< 90%) 
0.6833 (95-97.5%) 


1.2127 (85-90%) 
0.7817 (75-80%) 
2.4762 (97.5-99%) 



timescales (cf. §5.1.2) and the oldest stellar populations. Ob- 
viously, this is an oversimplified model, but it explains, at 
lest qualitatively, the observed trends. 



6 CONCLUSIONS 

This paper presents the study of a sample of 9 SO galaxies in 
the Fornax Cluster previously described in Paper I. By link- 
ing the properties of the central stellar populations (ages, 
chemical abundances) of these galaxies and their global 
properties (mass and dynamics), we uncover important clues 
on the main physical drivers governing the formation and 
evolution of SOs. Our main conclusions are the following: 

• Central absorption line indices correlate with central ve- 
locity dispersions in a way similar to what previous studies 
found for elliptical galaxies. However, a study of the stel- 
lar population properties of the Fornax SOs reveals that the 
trends shown by their line indices seem to be produced by 
relative differences in age and a-element abundances, con- 
trary to what is found in ellipticals where the overall metal- 
licities are important drivers of the correlations. 
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Figure 15. Central ages and metallicities from BC03 models and Mg6/(Fe) (oc o-elements overabundance) versus 3 different dynamical 
mass tracers. The same symbols are used as in Figure [14] 
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Figure 16. Central ages and metallicities from an improved version of Vaz96 models (Vazdekis et al. in preparation) versus 3 different 
dynamical mass tracers. The same symbols are used as in Figure [141 



• The scatter in the Index*-log(ao) relations can be 
partially explained by the rotationally-supported nature of 
many SOs. The tighter correlations found between Index* 
and log(l/MAx) support this interpretation. 

• The dynamical mass seems to be the primary physi- 



cal property governing these correlations and in the Fornax 
SOs we need to study it by considering their rotationally- 
supported nature. For these systems, Vmax is a better tracer 
of dynamical mass than a. 

• The a-element overabundance of these galaxies seems 
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Figure 17. Central ages and metallicities from BC03 models and 
Mgb/ (Fe) (oc a-elements overabundance) versus eachother. 

to be correlated with their mass. Moreover, the absorption- 
line-derived ages also correlate with the overabundances. 
This implies that the most massive SOs have the shortest 
star-formation timescales and the oldest stellar populations. 

These results support the idea of different star forma- 
tion histories for the central regions of bright and faint SO 
galaxies. Other authors have reached similar conclusions by 
studying carly-typc systems in the Coma Cluster (Poggianti 
et al. 2001; Mehlert et al. 2003). 

The results found here are consistent with a scenario 
where faint SOs arc descendants of spiral galaxies who lost 
(or exhausted) their gas and suffered a final central gasp 
of star formation during cluster infall, while bright systems 
seem to have star formation histories resembling those of 
normal ellipticals. However, are the nuclei characterising the 
entire galaxies? Have the outermost regions another story to 
tell? These questions lie beyond the scope of this study, but 
will be addressed in the following paper of the series, when 
a stellar population study at larger galactocentric distances 
will give us new clues about the detailed star formation his- 
tory of SOs in the Fornax Cluster. 
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APPENDIX A: TABLES 

In this appendix we include tables with different parameters 
for each galaxy of our sample. 

Table Al includes central absorption line indices mear- 

surod at 3 A resolution for all our galaxies while Table A2 
presents similar values measured at Lick resolution. 

In Table A3, central luminosity-weighted ages and 
[Fe/H] are presented for our galaxies by using different 
[Fe/H] indicators, H/3 as the age indicator and Bruzual & 
Chariot (2003) models at Lick resolution. 



Table Al. Central Line Indices measured at 3 A resolution of SO galaxies in Fornax. 



Name 


[MgFe]' 


(Fe) 


H/3 


Fe5015 


Mgi 


Mg2 


Mgb 


Fe5270 


Fe5335 


Fe5406 


Fe5709 


Fe5782 




[A] 


[A] 


[A] 


[A] 


[mag] 


[mag] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


CENTRAL 


NGC 1380 


4.14(0.04) 


3.53(0.03) 


1.48 (0.04) 


6.10(0.08) 


0.149(0.001) 


0.305(0.001) 


4.92(0.03) 


3.43 (0.04) 


3.64 (0.04) 


2.29(0.03) 


1.11 (0.03) 


1.02(0.03) 


NGC 1381 


3.71 (0.04) 


3.09(0.03) 


1.65 (0.04) 


5.46 (0.08) 


0.113(0.001) 


0.263(0.001) 


4.42 (0.03) 


3.14(0.04) 


3.04 (0.04) 


2.03 (0.03) 


1.16(0.03) 


0.87(0.03) 


NGC 1380A 


3.18 (0.09) 


2.99 (0.06) 


2.91 (0.08) 


5.68 (0.16) 


0.058 (0.001) 


0.183 (0.002) 


3.33 (0.08) 


3.10(0.08) 


2.88 (0.10) 


1.77 (0.07) 


0.98 (0.07) 


0.86 (0.07) 


NGC 1375 


2.71 (0.08) 


3.06 (0.06) 


3.89 (0.09) 


6.76 (0.17) 


0.053 (0.001) 


0.161 (0.002) 


2.43 (0.08) 


2.98 (0.08) 


3.15(0.09) 


1.78 (0.06) 


1.10 (0.05) 


0.76 (0.05) 


IC 1963 


3.54 (0.08) 


3.31 (0.05) 


2.45 (0.07) 


6.40 (0.14) 


0.072 (0.001) 


0.212 (0.001) 


3.76 (0.07) 


3.38 (0.07) 


3.24 (0.08) 


1.97 (0.06) 


1.16 (0.05) 


0.79 (0.05) 


ESO358-G006 


2.59 (0.16) 


2.61(0.11) 


2.35 (0.15) 


4.75 (0.29) 


0.039 (0.003) 


0.142 (0.003) 


2.51(0.14) 


2.75(0.15) 


2.48 (0.17) 


1.52 (0.13) 


1.01(0.11) 


0.58 (0.12) 


ESO358-G069 


2.82 (0.07) 


2.56 (0.05) 


2.25 (0.06) 


4.71 (0.12) 


0.048 (0.001) 


0.155 (0.001) 


3.03 (0.06) 


2.70 (0.06) 


2.42 (0.07) 


1.69 (0.05) 


0.94(0.05) 


0.46 (0.05) 


NGC 1316 


3.73 (0.08) 


3.28 (0.05) 


2.16(0.07) 


5.96 (0.14) 


0.094 (0.014) 


0.228 (0.001) 


4.26 (0.07) 


3.26 (0.07) 


3.30 (0.08) 


2.07 (0.06) 


1.04(0.05) 


1.02 (0.05) 


ESO 359-G002 


1.95 (0.15) 


1.92 (0.11) 


3.54 (0.13) 


3.99 (0.26) 


0.025(0.011) 


0.107 (0.003) 


2.01 (0.13) 


1.87 (0.14) 


1.97(0.16) 


1.12(0.12) 


0.57(0.11) 


0.63(0.11) 



Notes: From (2) to (13), la rms errors between () : Col (1), name; Col (2), [MgFe]' combined index (Gonzalez 1993; Thomas, Maraston & Bender 2003); Col (3), (Fe) combined index 
(Gorgas, Efstathiou & Aragon-Salamanca 1990); Col (4), H/3 index; Col (5), Fe5015 index; Col (6), Mgi index in magnitudes; Col (7), Mg2 index in magnitudes; Col (8), Mgb index; 
Col (9), Fe5270 index; Col (10), Fe5335 index; Col (11), Fe5406 index; Col (12), Fe5709 index; Col (13), Fe5782 index. 



Table A2. Central Line Indices measured at Lick resolution of SO galaxies in Fornax. 



Name 


[MgFe]' 


(Fe) 


H/3 


Fe5015 


Mgi 


Mg2 


Mg6 


Fe5270 


Fe5335 


Fe5406 


Fe5709 


Fe5782 




[A] 


[A] 


[A] 


[A] 


[mag] 


[mag] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


CENTRAL 


NGC 1380 


3.82(0.04) 


3.01 (0.03) 


1.47 (0.04) 


5.13(0.08) 


0.162'' (0.007) 


0.322" (0.007) 


4.86 (0.04) 


3.02 (0.04) 


3.00 (0.04) 


1.97 (0.03) 


1.00 (0.03) 


0.83(0.03) 


NGC 1381 


3.40 (0.04) 


2.70(0.03) 


1.65(0.04) 


4.83 (0.08) 


0.128" (0.007) 


0.283" (0.007) 


4.20 (0.04) 


2.81 (0.04) 


2.60 (0.05) 


1.75 (0.03) 


1.05 (0.03) 


0.75(0.03) 


NGC 1380A 


2.88 (0.09) 


2.65 (0.06) 


2.78 (0.08) 


4.69 (0.16) 


0.075" (0.011) 


0.204" (0.007) 


3.07 (0.08) 


2.79(0.09) 


2.51 (0.10) 


1.53 (0.07) 


0.90 (0.07) 


0.66(0.07) 


NGC 1375 


2.52 (0.07) 


2.63 (0.05) 


3.47 (0.07) 


5.24 (0.14) 


0.057" (0.011) 


0.169" (0.007) 


2.39 (0.07) 


2.69 (0.07) 


2.57(0.08) 


1.55 (0.00.1626) 


1.00 (0.05) 


0.68 (0.05) 


IC 1963 


3.22 (0.08) 


2.87 (0.05) 


2.26 (0.07) 


5.27 (0.14) 


0.089" (0.015) 


0.233" (0.007) 


3.56 (0.07) 


3.00 (0.07) 


2.75 (0.08) 


1.71 (0.06) 


1.02 (0.05) 


0.64 (0.05) 


ESO358-G006 


2.41 (0.16) 


2.23 (0.12) 


2.11(0.15) 


3.97 (0.29) 


0.054" (0.028) 


0.161" (0.008) 


2.51(0.14) 


2.44(0.15) 


2.03 (0.18) 


1.35 (0.13) 


0.89 (0.12) 


0.50 (0.12) 


ESQ 358-G069 


2.53 (0.07) 


2.25 (0.05) 


2.16(0.06) 


3.92 (0.13) 


0.065" (0.008) 


0.178" (0.007) 


2.76 (0.06) 


2.40 (0.06) 


2.10 (0.07) 


1.47 (0.05) 


0.87 (0.05) 


0.34 (0.05) 


NGC 1316 


3.22 (0.08) 


2.76 (0.05) 


2.11(0.07) 


5.15 (0.14) 


0.111" (0.015) 


0.250" (0.007) 


3.74(0.07) 


2.77 (0.07) 


2.75 (0.08) 


1.66 (0.06) 


0.93 (0.05) 


0.77 (0.05) 


ESO 359-G002 


1.70 (0.14) 


1.54 (0.10) 


3.32 (0.13) 


3.15 (0.25) 


0.039" (0.012) 


0.123" (0.008) 


1.86 (0.13) 


1.58 (0.14) 


1.51 (0.16) 


1.00 (0.12) 


0.52 (0.11) 


0.51 (0.11) 



Notes: From (2) to (13), Icr rms errors between () ; Col (1), name; Col (2), [MgFe] combined index (Gonzalez 1993; Thomas, Maraston & Bender 2003); Col (3), (Fe) combined index 
(Gorgas, Efstatliiou & Aragon-Salamanca 1990); Col (4), H/3 index; Col (5), Fe5015 index; Col (6), Mgi index in magnitudes; Col (7), Mg2 index in magnitudes; Col (8), Mgfe index; 
Col (9), Fe5270 index; Col (10), Fe5335 index; Col (11), Fe5406 index; Col (12), Fe5709 index; Col (13), Fe5782 index. 
(") An additive spectrophotometric correction to the Lick system was applied to the indices Mgi and Mg2. 



Table A3. Central Ages and Metallicities using line indices at Lick resolution for SO galaxies in Fornax. 



Name log(Age)[MgFe]' [Fe/HJjMgFe]' log(Age)^pe) [Fe/HJ^Fe) log(Age)Mgb [Fe/H]Mgb log(Age)FeB709 [Fe/HjFeSTOQ 

[yr] [dex] [yr] [dex] [yr] [dex] [yr] [dex] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



NGC 1380 


10.09 


+0.06 
-0.05 


0.23 


+0.05 
-0.05 


10.21 


+0.06 
-0.06 


-0.02 


+0.03 
-0.03 


9. 


.95 


+0.05 
-0.05 


0.63 


+0.06 
-0.06 


10.23 


+0.06 
-0.05 


-0.05 


+0.05 

-0.05 


NGC 1381 


Q 

y. 




+0.05 
-0.05 


0.06 


+0.06 
-0.05 


10.10 


+0.06 
-0.06 


-0.17 


+0.03 
-0.04 


9. 


.89 


+0.06 
-0.06 


0.34 


+0.07 
-0.06 


y. 


QQ 

.yy 


+0.04 
-0.04 


0.06 


+0.06 
-0.06 


NGC 1380A 


9. 


.16 


+0.06 
-0.06 


0.53 


+0.12 
-0.11 


9.23 


+0.04 
-0.04 


0.33 


+0.06 
-0.06 


9. 


.08 


+0.08 
-0.08 


0.77 


+0.16 
-0.15 


9. 


.30 


+0.04 
-0.04 


0.07 


+0.15 
-0.17 


NGC 1375 


8. 


.87 


+0.06 
-0.05 


0.67 


+0.10 
-0.10 


8.91 


+0.04 
-0.04 


0.56 


+0.04 
-0.04 


8. 


.84 


+0.06 
-0.07 


0.78 


+0.14 
-0.13 


8. 


.97 


+0.04 
-0.06 


0.36 


+0.08 
-0.09 


IC 1963 


9. 


.41 


+0.06 
-0.06 


0.42 


+0.10 
-0.10 


9.43 


+0.06 
-0.05 


0.27 


+0.06 
-0.05 


9. 


.36 


+0.07 
-0.07 


0.59 


+0.14 
-0.13 


9. 


.47 


+0.07 
-0.06 


0.14 


+0.11 
-0.11 


ESO358-G006 


9, 


.81 


+0.14 
-0.13 


-0.37 


+0.16 
-0.15 


9.79 


+0.13 
-0.12 


-0.30 


+0.12 
-0.12 


9, 


.83 


+0.14 
-0.13 


-0.43 


+0.14 
-0.13 


9. 


.72 


+ 0.15 
-0.14 


-0.10 


+0.25 
-0.30 


ESO358-G059 


9, 


.74 


+0.06 
-0.07 


-0.25 


+0.07 
-0.07 


9.75 


+0.06 
-0.06 


-0.26 


+0.05 
-0.05 


9, 


.74 


+0.07 
-0.08 


-0.23 


+0.07 
-0.06 


9. 


.70 


+0.06 
-0.06 


-0.13 


+0.11 

-0.12 


NGC 1316 


9, 


.54 


+0.09 
-0.10 


0.28 


+0.09 
-0.09 


9.61 


+0.07 
-0.08 


0.12 


+0.05 
-0.05 


9. 


.41 


+0.09 
-0.07 


0.57 


+0.13 
-0.13 


9. 


.69 


+0.08 
-0.08 


-0.04 


+0.11 
-0.12 


ESC) :!.')') c:i)i)L> 


9 


19 


+ 0.0,'') 

II. 0.1 


-0.27 


+ 0.16 

0.1."; 


9.2.3 


+ 0.05 

n.0.1 


-O.l.T 


+ 0.10 

11. iO 


9 


.13 


+ 0.06 

11.07 


-0.01 


+0.17 

il.i.1 


9, 


,37 


+ 0.17 
0.l:i 


-0.89 


+0.37 
-0.42 



Notes: From (2) to (9), 99% confidence intervals are presented. Col (1), name; Col (2) (3), age and metallicity, estimated using Bruzual & Chariot (2003) simple stellar population 
models at Lick resolution and [MgFe]' as metallicity indicator; Col (4) (5), age and metallicity, estimated using (Fe) as metallicity indicator; Col (6) (7), age and metallicity, estimated 
using Mg6 as metallicity indicator; Col (8) (9), age and metallicity, estimated using Fe5709 as metallicity indicator. 
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